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The previous treatments for strange quark matter in tiie quark mass-density-dependent model iiave 
unreasonable vacuum limits. We provide a method to obtain the quark mass parametrizations and 
give a self-consistent thermodynamic treatment which includes the MIT bag model as an extreme. In 
this treatment, strange quark matter in bulk still has the possibility of absolute stability. However, 
the lower density behavior of the sound velocity is opposite to previous findings. 
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I. INTRODUCTION 

Since Witten's conjecture that quark matter with 
strangeness per baryon of order unity might be bound 
[01, an extensive body of Hterature has investigated the 
stabihty and/or probabihties of strange quark matter 
(SQM) 1^. Because the apphcation of perturbative quan- 
tum chromodynamics (QCD) to strong-couphng domain 
is unbelievable while the lattice approach is presently 
limited to the case of zero chemical potential, we have 
to resort to phenomenological models. One of the most 
famous models is the MIT bag model with which Farhi 
and Jaffe find that SQM is absolutely stable around the 
normal nuclear density for a wide range of parameters 

. Further investigations have also been carried out by 
many other authors in the bag model |^,^,^ . A recent in- 
vestigation indicates a link of SQM to the study of quark 
condensates while a more recent work has carefully 
studied the relation between the charge and critical den- 
sity of SQM §. 

Chakrabarty et al. [^,|lO| have discussed the limita- 
tion of the conventional MIT bag model which assumes 
that the quarks are asymptotically free within the bag. 
In order to incorporate the strong interaction between 
quarks, one way is to fall back on the perturbation the- 
ory, which is questionable in the strong-coupling domain. 
An alternative way is to make the quark masses density- 
dependent. In this nonperturbative treatment, the strong 
interaction between quarks is mimicked by the proper 
variation of quark masses with density. There are two 
questions of crucial importance to this model. One is how 
to parametrize quark masses, the other concerns thermo- 
dynamic treatment. However, the two aspects are not 
self-consistent in literature presently. 

Here are the popularly used parametrizations for quark 
masses rriq [q — u,d, s,): 



mu,d = 



B 



B 



(1) 
(2) 



where m^o is the s quark current mass, rib is the baryon 
number density, B is the famous MIT bag constant. 
Equation (|^) was first used to study light quark matter 
| pT[ , and later extended to Eq. (^) to investigate strange 
quark matter ||,|o|jl|. 

As for the thermodynamic treatment, there exist two 
controversial ones in literature up to now. One expresses 
the total pressure of SQM as |||,|lO| 



Pi 



(3) 



where D, is the ordinary thermodynamic potential density 
of SQM [see Eq. @]. The other adopts the following 
expression ]l2[: 



P2 



drib' 



(4) 



The extra term in Eq. is said to arise from the baryon 
density dependence of quark masses. This difference 
leads to significantly different results. Therefore, it is 
meaningful to take a check of the two thermodynamic 
treatments. 

As is well known, the QCD vacuum is not necessarily 
empty. To obtain the vacuum properties, let us take the 
limit Tifc — !■ for the two treatments. It is easy to obtain, 
at zero temperature, the limits 



lim Pi 

nt—fO 



0, 



lim El = B, 

for the first treatment, and the limits 
lim Pa ^ -B, 
lim Ps = 2P, 



(5) 
(6) 

(7) 
(8) 



for the second treatment. Here Pi and P2 are the corre- 
sponding energy densities. 

According to the fundamental idea of MIT bag model, 
QCD vacuum has a constant energy density P, the fa- 
mous bag constant. The mass parametrization Eq. (m) 
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is just obtained from this requirement (lim„j_>o Ei — > 
SmqUb for flavor-symmetric case) The constant vac- 
uum energy comes from the fact that QCD vacuum must 
have a pressure to maintain pressure balance at the bag 
boundary. Obviously, the first treatment can give the 
correct vacuum energy and a wrong QCD vacuum pres- 
sure. On the contrary, the second treatment leads to 
the correct QCD vacuum pressure but a wrong vacuum 
energy. In fact, this is just caused by the ignorance of 
the QCD vacuum energy which guarantees the pressure 
balance at the bag boundary. 

It should be pointed out that in getting the unrea- 
sonable limits (||)-(||), we have used the quark mass 
formulas (0) and (||). Because these formulas are pure 
parametrizations without any real support from under- 
lying theories, one may ask if the contradictions can be 
solved by choosing other parametrizations? According to 
our present investigation, one should modify the quark 
mass formulas and thermodynamic treatment simultane- 
ously. 

It is the aim of this paper to give a self-consistent treat- 
ment which includes the conventional MIT bag model as 
an extreme. In our new treatment, strange quark matter 
in bulk still has the possibility of absolute stability. How- 
ever, the lower density behavior of the sound velocity in 
SQM is opposite to previous findings. 

In the following section, we first derive the new quark 
mass formulas and describe our thermodynamic treat- 
ment, and then in Sec. Ill, we present our results in 
studying SQM with this model. Section [V is a short 
summary. 



II. FRAMEWORK 

Let us schematically write the QCD Hamiltonian den- 
sity as 



QCD 



H, 



T-qoqq 



(9) 



where is the kinetic term, rrigo is the quark current 
mass, and Hi is the interaction part. The summation 
goes over all flavors considered. 

The basic idea of the quark mass-density-dependent 
model of strange quark matter is that the system en- 
ergy can be expressed as the same form with a proper 
noninteracting system. The strong interaction between 
quarks is included within the appropriate variation of 
quark masses with density. In order not to confuse with 
other mass concepts, we refer such a density-dependent 
mass to an equivalent mass in this paper. Therefore, if 
we use the equivalent mass m^, the system Hamiltonian 
density should be replaced by an Hamiltonian density of 
the form 



q 



(10) 



where vriq is the equivalent mass to be determined. Ob- 
viously, we must require that the two Hamiltonian den- 
sities iJoqv and -ffqcD have the same eigenenergy for any 
eigenstate l^"), i.e.. 



(*|iJeqv|*) = (*|i?QCD|*). 



(11) 



Applying this equality respectively to the state jrib) with 
baryon mmrber density rif, and the vacuum state |0), and 
then taking the difference, one has 



\nb) 



(0|i?QCD|0). 

(12) 



The simplest and most symmetric solution for the equiv- 
alent mass from this equation is 



= m,o + mi, 
where we have used the symbol definitions: 



(13) 
(14) 



{nb\Hi\7n,), {H, 



|i/i|0), and (q-g)„^ 



{nb\qq\nb) 



(qq), = {0\qq\0). 

Therefore, if quarks are decoupled, they should take 
the equivalent mass of the form (|3|) to keep the system 
energy unchanged. From Eq. ( p^ ) we see that the equiva- 
lent mass rrig includes two parts: one is the original mass 
or current mass mqo, the other is the interacting part 
nil. Because mi equals to the ratio of the total inter- 
acting part of the energy density and the total relative 
quark condensate, it is flavor-independent and density- 
dependent. Because of the quark conflnement and the 
asymptotic freedom, i.e.. 



lim mi 



oo, 



lim mi — 0, 

rib — 

the reasonable form might be 

D 

mi = — . 



Accordingly, we have 



D 

7^' 



(15) 
(16) 



(17) 



(18) 



where Z? is a free parameter to be determined by stabil- 
ity arguments. Obviously, z > for confined particles 
and z < for nonconfined particles. In Eqs. (|^) and 
(^), z = 1. However, just as mentioned in the Introduc- 
tion section, Eqs. (^ and (^ is closely linked to the first 
thermodynamic treatment, and thus unsuitable for our 
case. We now discuss the determination of mi which is 
consistent with our thermodynamic treatment. 
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Firstly, we express the interacting part of the energy 
density (Hi) = {Hi)n^ — {Hi)o [the numerator in Eq. (13)] 
as 

(Hi) = ^ J v{r){3nbd n)i3nbd r^) (19) 



(20) 



where r = | ri — r2 |, v(r) is the quark-quark interaction, 
R is the SQM radius, V = A/SttR^ is the volume. The 
extra factor 1/2 is responsible for double counting. 
Because of the following obvious equality: 



hm 

"6^0 {qq}Q 



1, 



(21) 



the Taylor series of the relative condensate at zero density 
has the following general form: 



1 



rib 



higher orders in nt 



(22) 



(99>o Pq 
If taking it only to first order approximation, we have 

E [(m)n, - (99)0] - E [-(^l)o/p'g] nb = Aub. (23) 



Taking the ratio of Eqs. ( |20[ ) and (^31), we get 



m\ = — T-nb \ v[r)r dr. 
A Jo 



(24) 



According to the lattice calculation |13| and string 
model investigation jl^], the quark-quark interaction is 
proportional to the distance, i.e., v(r) — ar. We thus 
have 



187ra i?^ 1 
mi = — - — rib—— oc 



(25) 



Therefore, we should take in Eq. ( p^ ) z = 1/3, i.e.. 



D 

iq — rriqQ -\ YJ^, 

nj 



(26) 



is the corresponding Fermi momentum. 

As usually done, we here assume the SQM to consist 
of u, d, and s quarks, and electrons (neutrinos enter and 
leave the system freely) . The degeneracy factor gi is 6 for 
quarks and 2 for electrons. The electron mass rUe is equal 
to 0.511 MeV. In order to include the strong interaction 
between quarks, the quark masses mq{q — u, d, s) should 
be replaced with the expression ([l3| ) or (p6|). The extra 
term B comes from the pressure balance condition, and 
its physical meaning is still the vacuum energy density 
or vacuum pressure just as in the MIT bag model. The 
corresponding pressure is 



P = 



^iHi - E, 



(29) 



i—u,d,s,e 



where is the chemical potential for particle type i. Be- 
cause it is equal to the Fermi energy at zero temperature, 
we have 



Pi 



Equation ( p9| ) is equivalent to 

p = -n-B. 



(30) 



(31) 



The second term — B is responsible for pressure balance. 
Such an extra term is necessary even in the nonrelativistic 
treatment of SQM ||l|l. 

It is clear that the above thermodynamic treatment 
will approach the conventional MIT bag model if one 
casts away the interacting part mi of the equivalent mass 



ruq. It can be proved, from Eqs. (p6[), (g7|), and (g9|), that 
we have the following correct vacuum limits: 

lim E = B, (32) 

lim P = -B. (33) 

Therefore, the physical meaning of B is the same as that 
in the conventional bag model. We take i?^/^ = 144 MeV 
in our present calculation. 



III. PROPERTIES OF STRANGE QUARK 
MATTER 



where Z? is a parameter to be determined by stability 
arguments. 

Because the Hamiltonian density -ffcqv has the same 
form as that of a system of free particles with equivalent 
mass rriq, the energy density of SQM can be expressed as 



where 



Pf.i 



p2 + mf p^dp + B, (27) 



1/3 



(28) 



Following previous authors 0, we assume the SQM to 
be a Fermi gas mixture of u, d, s quarks and electrons 
with chemical equilibrium maintained by the weak inter- 
actions: d,s'^u-\-e-\-Ve, s + u u -\- d. For a given 
baryon number density rib and total electric charge den- 
sity Q, the chemical potentials /i^, Ps, and /ig are 
determined by the following equations |8|: 

(34) 
(35) 
(36) 
(37) 



fid = Ps ^ P, 
Pu + Pe = P, 

i(7i„ + nd + ris) = rib, 
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where the particle number density is related to the 
corresponding chemical potential /i^ by 



9i I 2 2\3/2 



which is derived from the relation 

din 

with 



(38) 



(39) 



O. — Si 



3mfln ^-+^^ 



(40) 



In order to include the strong interaction between 
quarks, the quark masses m„,md, and nis in the above 
equations are to be replaced with the density-dependent 
expression Eq. (^6|) while the electron mass is negli- 
gible (0.511 MeV). 

For the bulk SQM with weak equilibrium, the previous 
investigations got a slightly positive charge. Our recent 
study demonstrates that negative charges could lower the 
critical density. However, too much negative charge can 
make it impossible to maintain flavor equilibrium. There- 
fore, the charge of SQM is not allowed to shift too far 
away from zero at both positive and negative directions. 
For this and our methodological purpose, we only con- 
sider neutral SQM in this paper, i.e., Q = in Eq. (^. 

Since the baryonic matter is known to exist in the 
hadronic phase, we must require D to be such that the ud 
system is unbound. This constrains D to be bigger than 
(47 MeV)2, i.e., at P = 0, Ejrn, > 930 in order not to 
contradict standard nuclear physics. On the other hand, 
we are interested in the possibility that SQM might be 
absolutely stable, i.e., at P = 0, E/ni, < 930, which gives 
an upper bound (128 MeV)^. we take D^/^ to be 50, 80, 
and 110 MeV, respectively. 

Because the light quark current masses are very small, 
their value uncertainties are not important. So we take 
the fixed central values rriuo — 5 MeV and m^o = 10 
MeV in our calculation. As for s quarks, we take 150, 
120, and 90 MeV, corresponding respectively to D^^^ = 
50, 80, and 110 MeV. 

For a given n^, we first solve for fii {i = u, d, s, e) from 
the equation group (|34|)~(|37|), and then calculate the en- 
ergy density and pressure of SQM from Eqs. (^7|) and 

E = j:^^^Fix,) + B^j:m,n,F{x,) + B, (41) 

i i 

P^E ^#^G(a;,) - 5 = E m.n.xf G(a;,) - B, (42) 

i i 

where the summation goes over u, d, s, and e, and 



rrii 



mi 



(43) 



is the ratio of the Fermi momentum to the mass that 
related to particle type i. With the hyperbolic sine func- 
tion sh^^(x) = ln(a; + Va;^ + 1), the functions F{x) and 
G{x) are defined as 

F{x) EE I [xVx^ + 1(2x2 ^ ^) _ sh-^a;)] /x^, (44) 
G{x) EE i [a;VPTT(2a;2 _ 3) + 3sh~^(x)] /x^, (45) 

which have the limit properties 

lim F{x) = 1, (46) 

x—*0 

limG(x) = i. (47) 



Therefore, we have the correct limits Eqs. (32) and (|32^ 

In Fig. ^ we give the energy per baryon vs baryon 
number density for the three pairs of parameters. We see 
that SQM is absolutely stable for the first two parameter 
groups, while metastable for the third group. The points 
marked with a circle "0" ar6 the zero pressure points 
where the pressure within SQM is zero. Because of the 
density dependence of quark masses, the zero pressure 
density is generally not that corresponding to the mini- 
mum energy per baryon (as in the usual case) , but nearly 
the case in the first two parameter groups. 

The resulting equation of state is plotted in Fig. ||. Be- 
cause it is insensitive to parameters, we have only chosen 
one parameter pair: D = (80 MeV)^ and niso ~ 120 
MeV. 

In Fig. ^, we show the sound velocity c of SQM with a 
dot-dashed line, which is obtained from 



dP 



dE 



1/2 



(48) 



Because the interacting part of the quark masses is neg- 
ligible at higher densities, it asymptotically tends to the 
ultrarelativistic value l/v^ as in the bag model (solid 
line). Simultaneously given is that calculated by the 
same method as in Ref. [||| with C = 90 MeV fm'^ 
and nisQ = 80 MeV (dotted line). Obviously, the lower 
density behavior of the sound velocity in our model is 
opposite to that in the previous calculation. 

It is interesting to note that if one considers the ther- 
modynamic relation P = —d{^lV) / dV as being more fun- 
damental than P = (as done in Ref. |l2)), Eqs. (|4^) 
and (^) should be replaced with 

E^J2 m,n,F{x,) + J2 m,nj{x^) + B, (49) 

i i 

P minixlG{xi) - min^f{xi) - B, (50) 



where 



fix.) 



3 rifc dnii 



2 TOj dub 
which has the limit property 



-h 1 - sh ^{xi) 



1^1 (51) 
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lim f{xi) = z. (52) 

However, the modification does not change the proper- 
ties of SQM significantly this time. For the same param- 
eters, the hne in Fig. 1 wiU move upward shghtly while in 
Fig. 2 and Fig. 3 a little downward. This is because the 
contribution from the extra term '^^minif{xi) (arising 
from the density dependence of the quark masses) is pos- 
itive to energy and negative to pressure. But the gross 
features of SQM are still the same. 

IV. SUMMARY 
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We have presented a new version of the quark mass- 
density-dependent model for SQM. We first note that 
the previous treatments have unreasonable vacuum lim- 
its. Then we provide a practical method to derive the 
quark mass formulas. In our thermodynamic treatment, 
the conventional bag model is included as an extreme, 
and the vacuum still has a constant energy density cor- 
responding to a constant pressure B. In this new treat- 
ment, SQM also has the possibility of absolute stability 
for a wide range of parameters. A noticeable feature 
is that the sound velocity is smaller than the ultrarela- 
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the points marked with circle "0-" 

FIG. 2. The equation of state for parameter group D^^^ = 
80 MeV and ruso — 120 MeV. It asymptotically approaches 
to the ultrarelativistic case as expected. 

FIG. 3. The sound velocity vs energy density. The 
dot-dashed line is calculated with the method in this paper, 
while the dotted line is calculated with the same method in 
Ref. Their lower density behavior is obviously opposite. 
The full line is the ultrarelativistic case. 
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